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Abstract 


A review of different approaches to describe anisotropic fractal media is proposed. 
In this paper differentiation and integration non-integer dimensional and multi-fractional 
spaces are considered as tools to describe anisotropic fractal materials and media. We 
suggest a generalization of vector calculus for non-integer dimensional space by using a 
product measure method. The product of fractional and non-integer dimensional spaces 
allows us to take into account the anisotropy of the fractal media in the framework of 
continuum models. The integration over non-integer-dimensional spaces is considered. In 
this paper differential operators of first and second orders for fractional space and non¬ 
integer dimensional space are suggested. The differential operators are defined as inverse 
operations to integration in spaces with non-integer dimensions. Non-integer dimensional 
space that is product of spaces with different dimensions allows us to give continuum 
models for anisotropic type of the media. The Poisson’s equation for fractal medium, the 
Euler-Bernoulli fractal beam, and the Timoshenko beam equations for fractal material 
are considered as examples of application of suggested generalization of vector calculus 
for anisotropic fractal materials and media. 

PACS: 64.60.Ak; 45.10.Hj; 62.20.Dc 


1 Introduction 

Fractals are measurable metric sets with non-integer dimensions (DEI- The basic property 
of the fractal is non-integer Hausdorff dimension that should be observed at all scales. The 
definition of the Hausdorff dimension requires the diameter of the covering sets to vanish. The 
fractal structure of real materials cannot be observed on all scales. This structure exists only 
for scales R > Ro, where R(> is the characteristic size of atoms or molecules of fractal media. 
Isotropic fractal materials can be characterized by the relation between the mass Md(Wb) of 
a ball region Wb of fractal medium, and the radius R of this ball in the form 



( 1 ) 
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The parameter D is called the mass dimension of fractal medium. The parameter I), does 
not depend on the shape of the region Wp, or on whether the packing of sphere of radius Ro 
is close packing, a random packing or a porous packing with a uniform distribution of holes. 
Anisotropic fractal materials can be characterized by the power-law relation for the mass of the 
parallelepiped region Wp in the form 

M d (Wp) = M 0 (j^ 1 (^j 2 \ min{L x , L y , L z } » Rq, (2) 

where the parameter cp, is non-integer dimension along A^-axis, k = 1,2,3, and L x , L y , L z 
represent three edges that meet at one vertex. The parameter oik describes how to increase the 
medium mass in the case of increasing the size of the parallelepiped along one axis, when the 
parallelepiped sizes along other axes do not change. The sum D — ai + a 2 + a 3 is called the 
fractal mass dimension of the anisotropic fractal medium. 

Using (JTJ) and (}2]) . we can define a fractal material as a medium with non-integer mass 
dimension. The non-integer dimension does not reflect completely all specific properties of 
the fractal media, but it is an main characteristic of fractal media and materials. For this 
reason, we assume that continuum models with non-integer dimensional spaces allow us to 
derive important conclusions about the behavior of the fractal media. 

The main ways of describing of fractal media can be conventionally divided into the following 
five approaches. 

1) Analysis on fractal approach: The first approach is based on the use of methods of 
” Analysis on fractals” pH0] E, 21 0 ;S] • Unfortunately a possibility of application of the ” Analysis 
on fractals” to solve differential equations on fractals |0] for real problems of fractal materials 
is very limited due to weak development of this area of mathematics to the present time. 

2) Fractional-differential continuum model approach: The second approach is based on the 
use of the fractional derivatives of non-integer orders with respect to space coordinates to 
describe some properties of fractal materials [S]. Therefore the correspondent models can be 
called the fractional-differential models. It has been suggested by Carpinteri and co-workers 
in m im 12] , where so-called local fractional derivatives are used, and then developed in 
Hams Ha ei na Hang. Unfortunately there are not enough differential equations with these 
fractional derivatives that are solved for various problems of fractal materials. It should be 
noted that the usual Leibniz rule does not hold [20] for derivatives of non-integer orders (and 
integer orders n f 1). It is a characteristic property of fractional derivatives. 

3) Fractional-integral continuum model approach. The third approach has been suggested in 
[21, 22j [23. :2U [25], 26, [2?] and it is based on application of continuum models of fractal media. 
These models can be called fractional-integral continuum model because the integrations of non¬ 
integer orders are used. The kernels of fractional integrals are defined by power-law density of 
states [27]. The orders of fractional integrals are equal to the mass (charge or other physical) 
dimensions of media. In these models, the density of states is applied in addition to the notion of 
distribution functions such as density of mass, density of charge. There are a lot of applications 
of these continuum models in different fields of mechanics and physics (see [2Zj and references 


2 


cited therein). These models have been applied by Ostoja-Starzewski in [28] [29, 1301 I3T1 132] . A 
generalization of fractional-integral continuum models for anisotropic fractal media has been 
suggested by Ostoja-Starzewski and co-workers in [33], 34] EH 36, [37] EB, 32] • In these models, 
the differential operators are modified by the density of states. However, these operators have 
integer-order differential operators. 

4) Fractional space approach: The fourth approach uses the concept of a fractional space, 
which is characterized by non-integer (fractional) powers of coordinates. The fractional space 
approach has been suggested in papers [40] ;4l] 02] an d then it is used for applications in dif¬ 
ferent areas [13], 04] 05] 06] [27] (see also [07] 08] 09] [50]). This approach has been developed 
by Calcagni in [51j [52, 03] [54] (see also [55] [56]), and then it was generalized for anisotropic 
case by using a multi-fractional space (’’multi-scale space”) in [57]ESE2]• The first interpreta¬ 
tion of the fractional phase space is connected with fractional dimension space. The fractional 
dimension interpretation follows from the formulas for dimensional regularizations and it was 
suggested in [40]. In the paper [41, 02] the second interpretation of the fractional phase space 
is considered. This interpretation follows from the fractional measure [40j of phase space that 
is used in the fractional integrals, i.e. the integrals of non-integer orders. In the third inter¬ 
pretation, the fractional phase space is considered as a phase space that is described by the 
fractional powers of coordinates and momenta. In addition almost all Hamiltonian systems 
with fractional phase space are non-Hamiltonian dissipative systems in the usual phase space. 
It allows us to have the fourth interpretation of the fractional phase space as a phase space of 
power-law type of non-Hamiltonian systems. Using fractional space approach we can consider 
wide class of non-Hamiltonian systems as generalized Hamiltonian systems. Differentiation in 
fractional space approach can be used in two forms: (a) the usual derivatives with respect to 
fractional powers of coordinates [ 4011411, 42]; (b) the fractional derivatives of non-integer or¬ 
ders (fractional derivatives) with respect to coordinates [51], 52] E3[, [54]. The term ’’fractional 
space” is sometimes used for non-integer-dimensional space. This leads to confusion and mis¬ 
understanding. We use the term ” fractional space” for effective space M” with coordinates that 
are non-integer powers of coordinates of physical space. In the fractional space approach, the 
integer-dimensional spaces M n , the integration and differentiation of integer-orders for these 
spaces are used. The coordinates of fractional space are considered as effective coordinates 
that are fractional powers of real space coordinates of physical system. Note that we also can 
use effective spaces for non-integer-dimensional physical spaces. Expressions for the effective 
coordinates for fractional and non-integer-dimensional spaces differ by factors in the density of 
states. 

5) Non-integer-dimensional space approach: The fifth approach is based on application of in¬ 
tegration and differentiation for non-integer-dimensional spaces. The integration in non-integer 
dimensional space is well developed [00, SB 112], and it has a wide application in quantum held 
theory. The axioms for integrals in D-dimensional space is suggested in [60]. This properties 
are natural and necessary in applications [62]. Integration in D-dimensional spaces with non¬ 
integer D is used for dimensional regularization in quantum field theory [63] 011, 62] and in 
physical kinetics [65] 06]. Dimensional regularization is a way to get infinities that occur when 
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one evaluates Feynman diagrams in quantum theory. Differentiation in non-integer dimensional 
space is considered in [61] [67, OB]- In the papers [SI, 67] it was offered only a scalar Lapla- 
cian for non-integer dimensional space. Unfortunately the gradient, divergence, curl operator 
and the vector Laplacian [104] are not considered in [61] [67]. The scalar Laplace operators, 
whcih are suggested by Stillinger in [6T] and Palmer, Stavrinou in [67] for non-integer dimen¬ 
sional spaces, have successfully been used for effective descriptions in different areas of physics 
and mechanics such as quantum mechanics (see m Eij, m - da, m - izsi, m - isej, 
[87] - [91] , the diffusion processes [92], the general relativity [93J [94J , and the electrodynamics 
[95 , 96l [971195119911100, T01J. All these applications are based only on two generalization of the 
scalar Laplacian that are suggested in [611, [57j. To expand the range of possible applications 
of models with non-integer dimensional spaces it is important to have generalization of differ¬ 
ential operators of first orders (gradient, divergence, curl operators) and the vector Laplacian. 
The continuation in dimension is recently suggested in [68] [69] to define the gradient, diver¬ 
gence, curl operator and the vector Laplacian for non-integer dimensional space. It allows us 
to describe isotropic fractal media in the framework of continuum models with non-integer di¬ 
mensional spaces. To generalize non-integer dimensional space approach for anisotropic fractal 
media we can suggest to use the product measure approach suggested in [4D, SH S2] and [67]. 
Generalizations of the gradient, divergence, curl operators and the vector Laplace operator for 
non-integer dimensional and fractional spaces to describe anisotropic fractal media are not con¬ 
sidered by the product measure approach in [61] [67] and other papers. These generalizations 
are suggested in this paper as an extension of approach proposed in [62] 69] . 

To present more clearly some of the differences between these five approaches to describe 
fractal media distributed in the space R n , we present the following table. 


Approach 

Set / Space 

Integration 

Differentiation 

Analysis on fractal 

Fractal set 

Integration for 
fractal set 

Differentiation for 
fractal set 

Fractional-differential 
continuum model 

Integer-dimensional 
space R n 

Integrals for R n 

Fractional-order 
derivatives for R n 

Fractional-integral 
continuum model 

Integer-dimensional 
space R n 

Fractional-order 
integrals for R n 

Integer-order 
derivatives for R n 

Fractional space 

Integer-dimensional 
(effective) space R” 

Integer-order 
integrals for R n 

Integer-order 
derivatives for R n 

Non-integer- 
dimensional space 

Non-integer¬ 
dimensional space 

Integrals for 
non-integer¬ 
dimensional space 

Derivatives for non- 

integer-dimensional 

space 


Let us note some advantages and disadvantages of the the third, fourth and fifth approaches. 
The forms of the functions that define the density of states in the third and fourth approaches 
have arbitrariness due to the existence of different types of fractional integrals. In the fifth 
approach, the form of density of states is uniquely fixed by the expression of volume of the 
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region in the non-integer-dimensional space. In the third approach practically all properties 
of fractal media are reduced to coordinate transformations and to space curvature of power 
type. The transition by transformations to an effective Euclidean space cannot be considered 
as a consistent approach to describe fractal media. These transformations can be considered 
only as a part of mathematical method to solve some equations in the framework of continuum 
models with non-integer-dimensional space, i.e., in the fifth approach. Connections between the 
effective coordinates and physical coordinates for the non-integer-dimensional space approach 
are uniquely defined. In addition, the fifth approach allows us to use fractional derivatives 
and integrals to describe nonlocal type of fractal media that cannot be described by other 
approaches. 

In this paper, we consider the two last approaches based on the fractional space and the 
non-integer dimensional space to describe anisotropic fractal materials and media. We sug¬ 
gest a generalization of vector calculus for non-integer dimensional space that is product of 
spaces with different dimensions. In section 2, we discuss the product measure method to 
describe anisotropic fractal media. In section 3, the integration for fractional and non-integer- 
dimensional spaces are considered. In section 4, differential operations of Erst and second 
orders for fractional space and non-integer-dimensional space are suggested. In section 5, we 
give some examples of application of suggested generalization of vector calculus for anisotropic 
fractal materials and media. 


2 Product measure method 

2.1 Product measure for fractional and non-integer-dimensional spaces 

The product measure method can be applied to the non-integer dimensional spaces and to the 
fractional spaces: 

(I) Product measure for the fractional spaces. The product measure approach for the frac¬ 
tional spaces has been suggested in [44], 42], where fractional phase space is considered with its 
interpretation as a non-integer (fractional) dimensional space. In papers [41], 02], the following 
measure is used for generalized coordinates and momenta 

dfia(x k ) = c(a) \x k \ a ~ l dx k , (3) 

where the numerical factor c(a) is 

c(a) = 1/T(a). (4) 

We use the factor c(a) in the form (J3J) to get a relation with the Riemann-Liouville fractional 
integrations of non-integer order a. In the papers [40] [41, 02], it has been shown that the 
integration for the suggested fractional space is directly connected with integration for non¬ 
integer dimensional space up to numerical factor. Therefore the fractional space has been 
interpreted as a non-integer dimensional space. The suggested fractional space approach has 
been used in 03! 00] 05] 06] 07] for configuration and phase spaces, and then it has been applied 
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by Calcagni in [51, [52] [53:, [53] for space-time. The differentiation and integration in fractional 
space are considered as differentiation and integration with respect to non-integer (fractional) 
powers of coordinates. The differential operator of the first order is defined |3EE HU 12] S3) 41 > IS] 
by 

n 9 1 9 

a ’ k dQ k a\x k \«~'dx k ' U 

where Q k = sgn(x) \x\ a . The product measure approach also used in [33:, EH .35, |36> 39] to 
describe fractal media, but the fractional space as space with non-integer powers of coordinates 
has not been considered. Instead of the Riemann-Liouville fractional integration, which is used 
in (201 HU HU HU HU HU HE]) the product measure approach is used in (33] EU 135] 136] [39] for 
so-called modified Riemann-Liouville integrations in the integer dimensional spaces. 


(II) Product measure for the non-integer dimensional spaces. The product measure approach 
has been suggested by Palmer and Stavrinou in ra for non-integer dimensional spaces, where 
each orthogonal coordinates has own dimension. The methods of this paper can be considered as 
a modification of Stillinger’s [El] and Svozils’s [103] methods for the product measure approach. 
In the paper |67j . the product of the following single-variable measures is used 

d/j, ak (x k ) = c(a k ) \x k \ ak ~ x dx k (6) 


with the numerical factor 


c(a ) = 


2 7T Q / 2 


T(«/2) 

Note that c(a) is equal to the surface area S a (R ) of a-sphere with the radius R — 1, where 

2n a/2 

S a (R ) = — — R a 


(7) 


T(«/2) 


( 8 ) 


The integration for non-integer dimensional space by product measure approach is described 
in [67]. The scalar Laplacian operator for non-integer dimensional spaces is also suggested in 
[67] in the form 


S ^(a) 


/ d 2 a k — 1 3 \ 

“ \dx 2 k Xk dx k ) 


(9) 


where ( a ) = (ai, a 2 > ^ 3 ) is tl ie multi-index. Unfortunately definitions of the gradient, diver¬ 
gence, curl operations, and the vector Laplacian are not considered in [67]. 


2.2 Type of anisotropic fractal media described by non-integer di¬ 
mensional spaces 

In this section, we describe a type of anisotropic fractal media that can be considered by 
approach based on non-integer dimensional spaces. 
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Let us defined the parameter (k = 1; 2; 3) that describes the scaling property along 
Xfc-axis by 

O 1 Ol x — D dyzi (10) 

where D is a non-integer mass dimension of the material. Here d yz is the non-integer dimension 
(for example, the box-counting dimension) of the Y Z-cross-section, which is perpendicular to 
the A"-axis. We should assume that non-integer dimension d yz of the Y Z-cross-sections is the 
same for all points along the X-axis ( d yz (x) = d yz = const). If this condition is not satisfied, 
then we have a variable order ot\{x) that depends on the coordinates. 

In general, the parameters d xy , d xz , d yz cannot be considered as dimensions of the boundaries 
of fractal media. The boundary of fractal medium can be fractal surface with dimension D s . In 
the general case, D s can be greater than two. The non-integer dimensions of the cross-section 
are always less than two, 


0 < dxy < 2, 0 < d xz < 2, 0 < d yz < 2. 


( 11 ) 


Similarly (ITUlh we can define other two parameters ck 2 = oiy and «3 = a z . 

Let L 1; L 2 , and L 3 be the basis vectors that dehne a three-dimensional parallelepiped in 
M 3 . The parallelepiped region is the convex hull for these vectors 

{ n=3 

a,k Lfc : 0 < a fc < 1 

k =i 

For simplification we will consider the rectangular parallelepiped only. Anisotropic fractal ma¬ 
terial can be characterized by the power-law relation for the mass Mp(Wp) of the parallelepiped 
region (fl2l) in the form 

where the parameter oik is non-integer dimension along X^-axis, k = 1,2,3, and Lp. = | L* | 
are the magnitudes of vectors L&, k — 1,2, 3. The values Lk can be considered as three edges 
that meet at one vertex, and Rq is a characteristic size of particles (atoms or molecules) of 
fractal medium. The parameter oik describes how to increase the mass medium in the case of 
increasing the size of the parallelepiped region along one axis, when the parallelepiped sizes 
along other axes do not change. The sum D = + a 2 + a 3 is called the dimension of the 

anisotropic fractal medium. 

In general, the parameters a^. > 0 (k = 1, 2, 3) can be either less than unity or greater than 
unity. In all cases, the following conditions should be satisfied 


( 12 ) 


0 < cki + ct 2 + «3 = D < 3. (14) 

For example, the conditions 0 < a*, < 1 for all k G {1,2,3} hold for fractal media similar to 
3D Cantor dust. We assume that the parameter cka, > 1 describes a fractal flow and motion 
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of medium in AVdirection. This possibility is based on an assumption that trajectories of the 
medium particles in the A^-direction are fractal curve with the dimension a k > 1. For example, 
the Koch curve has a k = ln(4)/ln(3) ~ 1.262. We also assume that a k > 1 can be used for 
materials consisting of fractal molecular curves or fractal chains 


2.3 Single-variable measure 

For the fractional space approach and the non-integer dimensional space approach, we can use 
the single-variable measure 

d/i(ct, x) = c(a) |x | a_1 dx, (15) 

where a > 0 is a parameter that will be considered as a non-integer (fractional) dimension of 
the line, and c(a) is a function of a. Here, we take the absolute value of x in Ixj " -1 to consider 
positive and negative values of x. For a = 1, the numerical factor c(a) must be equal to 1 in 
order to have 

d/i( 1, x ) = dx. 

Using the product measure approach for R 3 with point coordinates x\, X 2 , £ 3 , the single-variable 
measures are 

d/x(a k , x k ) = c{a k ) \x k \ ak ~ l dx k . (16) 

In the product measure method we can use the following two ways to define a numerical factor 
c(a k ). In the first way of description, the factor is defined by a connection with integrals of 
non-integer orders a k in the integer dimensional space. In the first way, the factor is defined 
by a connection with integrals in spaces with non-integer dimensions a k along the A*,-axis. 

Let us give the effective coordinates for these two cases. 

(1) The fractional space approach is based on the use of the following new (effective) coor¬ 
dinates 


Qk Qkis^ki-^k) 


sgn(x fe ) \x k \ ( 


F(«fc + 1) 

that is connected with the single-variable measure of the form 

1 


dji{a k , Xj^j dQ k 


r(a fc ) 


\x k \ ak 1 dx k , 


where the numerical factor in the density of states is 

1 


c(a) = 


r(<*) 


(17) 


(18) 


(19) 


This form of c(a k ) is based on the connection of the Riemann-Liouville integrals of non-integer 
orders a k . 

(2) For non-integer dimensional space approach, we can use the effective coordinates 


7 r 


oik /2 


2TK/2 + 1) 


k X k (ot k ,X k ') 


sgn(x fc ) \x k \ ak , 


( 20 ) 






that is connected with the single-variable measure [6T1167] of the form 





(21) 


Here we take the density of states in the form 



( 22 ) 


such that the area of the sphere with radius r is equal to 



(23) 


The absolute values \xk\ can be interpreted as radii rk = \xk\ of sphere with non-integer dimen¬ 


sion cx fc . The presence of a factor of 2 for S a _i in (123]) is due to the fact that for a = 1, the 
variable r is integrated from —R to R , and when the limits are taken as 0 and A, one gets a 
factor of 2. 

The space with coordinates (1201) and the product of single-variable measures (1211) can be 
considered as a non-integer dimensional space. The parameter D = a± + a 2 + «3 can be 
interpreted as a dimension of the space. For aq = a 2 = cx 3 = 1, we get D = 3, i.e., the 
dimension of space is the usual integer dimension. If cq = oi 2 — oq — ck, where 0 < a < 1, we 
have a non-integer dimensional space for isotropic fractal materials. Regardless of the isotropic 
or the anisotropic case, the dimension of the space is given by 


D = cki + a 2 + «3- 


We have a non-integer dimensional space if at least one of the parameters ak is not equal to 1. 

2.4 Density of states 

A connection between the single-variable measure d/j,((Xk,Xk) of non-integer dimensional space 
and the measure d/i(l,Xfc) of integer dimensional space is 


Xfc') C\ (cT/c, Xfd) dfl (1, Xjd) 


(24) 


The functions Xk) should be considered as a density of states of fractal material [27] along 

the Xfc-axis. We can define the density of states along the A^-axis by the equation 



(25) 


where c 3 (D, x, y, z) is the density of state in the volume of material, and C 2 (d xy , x, y ), c 2 (d xz , x, z ), 
C 2 (d yz ,y,z) are density of states of the XY, XZ , YZ- cross-sections respectively. 
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The interpretation of the functions cs(D, x, y, z),C 2 (d xy , x, y), C 2 (d xz ,x, z), C 2 (d yz , y, z), and 
c\(ak,Xk), (k = 1,2,3), as densities of states has been suggested in [2?]. The density of states 
c n describes how closely packed permitted states of particles in the space R n . In the book 
[27] the density of states is defined by integrations of non-integer orders that is also called the 
fractional integration. For Xk G [a*,;^] the function ci(a k ,x k ) is defined in [27] by 


ci (a k , x^) 


1 

fw 


| Xk &k 


(26) 


This form is connected with the Riemann-Liouville fractional integral. In the paper [33], an 
expression for Ci(ak,Xk) has been suggested in the form 


Cl {Oikj Xk) Oik \^k Xk 


(27) 


that does not contain the gamma function. This form is connected with the modified Riemann- 
Liouville fractional integral. 

In this paper, we use integration in non-integer dimensional space instead of fractional 
integration. Then we should use the density of states in the form 

7]-“fc/ 2 

Cl(ff ‘- Il) = rM) hl (28) 

that is defined by the measure for integration in non-integer dimensional space. An application 
of density of states in the form (1281) allows us to get the expression for lengths 

r R r R 27t "/ 2 W 2 

/ M a i x ) = 2 M a , x) = R a = R a = VJR) (29) 

J-r Jo aT(a/ 2) T(a/2 + l) 

that coincides with the well-known value for non-integer dimensional volume. 


3 Integration in fractional and non-integer dimensional 
spaces 

3.1 Product spaces and product measures 

The integral for non-integer dimensional space is defined for a single-variable in [103]. It is 
useful for integrating spherically symmetric functions only. We can consider multiple variables 
by using the product spaces and product measures [67]. 

Using a collection of n = 3 measurable sets (144, /W D) with k — 1,2, 3, we form a Cartesian 
product W = W\ x 144 x 144 of the sets 144. The definition of product measures and an 
application of Fubini’s theorem gives a measure for the product set W = W\ x 144 x 144 as 

71=3 

M4U) = (p-cn X Ra 2 X /ia 3 )(4U) = \\p{oi kl W k ). (30) 

k =1 
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Then integration over a function / on W is 


n =3 


f(x 1 ,x 2 ,x 3 )dn B = 


'w 


<Wi JW 2 JWi 


f(xI,x 2 ,x 3 ) d/x(a k ,x k ). 


k =1 


(31) 


In this form, the single-variable measure may be used for each coordinate x k , which has an 
associated non-integer dimension a k , by the equation 


dfj,(a k ,x k ) = Ci(a k , x k ) dx k , (.k = 1,2,3), 
where Ci(a k ,x k ) is the density of states of the form 

Note that we use ci(a k ,x k ) without the factor 2. 

Then the total dimension of W = W± x W 2 x W 3 is D = aq + a 2 + a 3 . 


(32) 


(33) 


3.2 Reproduce the single-variable integration 

Let us reproduce the result for the single-variable integration in the form 

r 2i t d/2 C°° 

J w x s) = y( D /2) J 0 ^ r dr ( 34 ) 

for spherically symmetric function f(xi,x 2 ,x 3 ) = f[r ) in W\ x W 2 x IT 3 , where r 2 = (xi) 2 + 
(x 2 ) 2 T (1C3) 2 . For this function, we can perform the integration in spherical coordinates (r, (ft , 6). 
The Cartesian coordinates (xi,x 2 ,x 3 ) can be expressed by the spherical coordinates (r, 99 , 6), 
where r G [0, 00 ), (p G [0, 27t), 6 G [0, 7 r], by: 


X\ — r sin 6 cos 99 , 

x 2 = r sin 0 sin c/ 5 , 
x 3 = r cos 6. 

In this case, equation (l3Tj) becomes 


d^B f(xi,x 2 ,x 3 ) = 


<w 


(35) 

(36) 

(37) 


= A(a 1 ,a 2 ,a 3 ) [ dx 1 [ dx 2 [ dx 3 \x 1 \ ai 1 \x 2 \ a2 l \x 3 \ a3 1 f(x u x 2 , x 3 ) = 
Jw 1 Jw 2 Jw 3 


— d(aq, a 2 , a 3 ) 




de J(r, 6) 


01 + 02 + 03—3 . 
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where 


cos^l" 1 - 1 1 sin<^|" 2_1 | sin 0|" 1+ " 2_1 | cos^ 3 " 1 /(r), 

r^Oi\j2 yj- (%2 /2 j^CX.3/2 

A (“^ " 2 ’ " s) = r(a 1 /2)r(a 2 /2)r(a 3 /2)’ (3§) 

and J(r, 0) = r 2 sin 0 is the Jacobian of the coordinate change. 

Since the function is only dependent on the radial variable and not the angular variables, 
we get the product of three integrals 


dfjiB f{xi,x 2 ,x 3 ) = A(a 1 ,a 2 ,a 3 ) / f(r) r “ 1+ " 2+Q3 1 • 


>w 

r 2n 


I dy\cosip\ ai ~ l Isin^l" 2 - 1 / dd | sin 0|" 1+Q2_1 | cos#!" 3-1 

' 0 Jo 


Using equation (26) from Section 2.5.12 of [105]. in the form 

r-K/2 


(sin xY - 1 (co sx)- l dx = f[!‘( 2) r( ( /2) . 

2T{(n + v)/2) 


where /1 > 0, u > 0, we have 


-2ir 


dtp\costp\ ai 1 |sin(^| 


CX-2 — 1 _ 


= 4 


t/2 


= 2 / d 6 (sin 0) 

Jo 


T/4 , , \cti —1 / • ^,-1 2 T(a 1 /2) r(a 2 /2) 

^ (C08 ^ ^ = r(( ai +a;)/2) ' 

/*7T 

/ d6 | sin#| ai+a2_1 1 cos0| a3_1 = 

a 3 —1 _ r ((«i + a 2)/2) r(a 3 /2) 


'0 


ai+Q2 —1 


(cos 6*) 


r((cKi + «2 + a 3 )/ 2 ) 

Using ZJ = O'! + a 2 + a 3 , we obtain for f(xi,X2,x 3 ) = f(r) the ralation 


(39) 


(40) 


(41) 


(42) 


r 2 7r D/2 Z 00 

/ f{xi,X2,x 3 )dfj, 1 (x 1 )dfj,2(x2)dfj, 3 (x 3 ) = r , n/0 c / 

1 1 -^ 7 U Jo 


f(r)r u 1 dr 


(43) 


This equation describes the .D-dimensional integration of a spherically symmetric function, and 
reproduces the result (IMP . 

ft is important to note that relation (|43|) holds only for density of states c(ak,Xk ) that 
corresponds to the non-integer dimensional space (l33jh Equation (1431) cannot hold only for 
density of states ci(ak,Xk) suggested in [22] and [33] . 
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4 


Differential operators for fractional and non-integer 
dimensional spaces 


4.1 Laplace operator for non-integer dimensional and fractional spaces 


Let us consider possible definitions of the scalar Laplace operators for non-integer dimensional 
and fractional spaces. 

(I) Scalar Laplace operators for non-integer dimensional space. The Laplacian operator for 
non-integer dimensional space has been suggested in [67] in the form 


S^(a) 


vf— + ak ~ 1 d A 

f^\dx 2 k x k dx k /’ 


(44) 


where (a) = (ai,a 2 ,a 3 ). Unfortunately definitions of differential operators of first order such 
as gradient, divergence, curl operations are not considered in [67]. The Laplacian operator (1441) 
is used in [98l [99 . lOOi 4011 [102] , In the book [102] . the gradient, divergence, and curl operators 
are suggested only as approximations of the square of the Laplace operator (J9]). The first order 
differential operator proposed in [402] is considered as an approximation of (144]) in the form 


D 


otkik 


d a k ~ 1 

2 x k 


dx 


For example, the gradient is 


grad D ip « 

k =1 


(dp a k - 1 
1 O- 1 -o-^ ) efc ’ 

V dx k 2 x k 


and the divergence for vector field u = u k e k is defined by 

,. f ® Uk 1 ak ~ 4 

d,VDU “U(a + 

k =1 


(45) 


(46) 


(47) 


Obviously, the corresponding scalar Laplacian 


div^grad D p 


3 

E 

k =1 


d 2 p Qi k ~ 
dxl x k 


1 dp (a k - l)(a x - 3) 


dxi 


+ 


4 4 


-p 


(48) 


does not coincide with the operator (1441) . 

(II) Scalar Laplace operators for fractional space. The scalar Laplace operators have been 
suggested by Calcagni 0 for for fractional space-time. We can represent the suggested equa¬ 
tions for fractional space in the following forms 


«iV>=E 


d 


—' v(a,x) dx k 


v(a, x) 


dp 

dx k 


(49) 
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V 2f = Y. 


a 2 


k =1 


and 


y/v(a,x) dx 2 k 

3 ( x k ) l ~ 1/2 <9 


(y/v{a,x) v), 


^ = g 7^3^ ( (I ‘ ),_1/2 ^ ( W '" I/2 WThy) 

where Xk > 0, and v(a,x) is the isotropic measure weight 


(50) 


(51) 


v(a,x) = JJci(a,Xfc) = Jj 


(a:*,) 


a— 1 


(x k > 0). 


(52) 


r(a) ’ 

k=1 k =1 v ' 

Here ci^a^Xk) is the density of states for fractional space [27]. We should note that Calcagni 
consider the Laplace operators for Minkowski space-time Mf 3 . Equations (1491) . ([50]) and (T5TT) 
are given for the Euclidean space M 3 to describe fractal material. Note that the expression (T5T1) 
contains (T49|) and (150]) as particular cases. The operators (1491) . (I50|) can be represented as 


Substitution of 


into 


/Cl — fc'a,l=l—a/2i /C2 — ^a,1=1/2- 

), (154]) gives the Laplace operators in the forms 
3 


JC*P = 


k=1 


d 2 p a — 1 dip 
dx\ x k dx k 


K.^ = /T 


k =1 
3 


d 2 p a — 1 dp (a — l)(a — 3) 
dx\ + x k dx k + 4 x\ 


f d 2 p a — 1 dp (a — 2) 2 — 4 1 


fc'ajtp 'y ] I n 2 4" rj 

\ dx k x k dx k 


+ 


4 4 


-p 


(53) 

(54) 

(55) 

(56) 


The Laplace operator (154]) coincides with the expression (jUJ) suggested in [67]. The Laplace 
operator (1551) coincides with the operator (1481) that can be derived from the first order operators 
(1451) and (145]) suggested in [ 102] . It was proved [54] that the case l = 1/2 is unique because it 
is only one, where the Laplace operator of the type (1541) can be represented as the square of 
first order differential operator. This Erst order operator (see equation (3.43) in [54] ) is 


T^a,kP 


d 


such that 


■\/v(a, x) dx k 
3 

IC 2 V = 

k= 1 


^V' u ( Q! > a; ) ¥>)> 

(57) 

\ 2 


,k J ' 

(58) 
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Substitution of (1521) into (j57li gives 


Tda,k^P 


dip a — 1 
dx k 2 x k 


(59) 


We can see that operator (IMJD can be considered as (1T5|) with a k = a for all k = 1,2, 3. 

To anisotropic fractal media and materials in the framework of the non-integer dimensional 
space approach we should generalize the isotropic measure weight v(a, x ) by using 


3 

v(at,x) = 

k=1 


7j- Q fc/ 2 

r(a fc /2 + 1) 


^r - 1 


(60) 


instead of the expression (|52l) for fractional space. The derivative of the first order (|59li should 
also be generalized for anisotropic fractal media as 


'Da,k l p 


dip 

dx k + 


Oik - 1 
2 x k 


<p, 


(61) 


where we take into account different values of dimensions along the X^-axis. Using the first 
order differential operator (1BTT) . it is easy to define the del operator, gradient, divergence, curl 
operators, and the vector Laplacian in order to describe fractal media in the framework of the 
non-integer dimensional space approach. 


4.2 Approaches to formulation of the vector calculus in non-integer 
dimensional space 

In the Stillinger’s paper [El], the scalar Laplace operator for non-integer dimensional space is 
suggested only. Generalizations of gradient, divergence, curl operators, and the vector Laplacian 
are not considered in [6T]. A generalization of the gradient, divergence, curl operators, the scalar 
and vector Laplace operators for non-integer dimensional space can be defined by the method of 
continuation in dimension [68l 169 |. For simplification, only radial case is consider in [68] . where 
the scalar and vector fields are independent of the angles, and the vector fields are directed 
along the radius vector. This simplification is analogous to consideration of the integration in 
non-integer dimensional space in Section 4 of [62] . 

The main advantage of the product measure approach to define the vector calculus for 
non-integer dimensional space is a possibility to describe anisotropic fractal materials. In this 
paper, we suggest differential vector operators for non-integer dimensional space by product 
measure method to describe anisotropic fractal media in the framework of continuum models. 
The differential operators are defined as an inverse operation to integration in non-integer 
dimensional space. 

We can state that it is possible to reduce non-integer dimensional space to a power-law 
curved space by considering the density of states Ci(ak,x k ) as the Lame coefficients 

H k = ci(a k ,x k ). (62) 
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We consider Euclidean space for the effective coordinates X k (x k ), which are defined by (l20]h 
such that 


d 2 s x = Y,(dX t f = yV|KH:.) (dx t ) 2 . 


(63) 


k= 1 


k=1 


where C\(a k ,x k ) are the density of states of the form 

7j-“fc/2 

Cl( "‘' I ‘ )= rW2 ) |ltl '‘" 1 

Using (163T) . we can see that the densities of states fl64|) are the Lame coefficients 


Hu = 


(64) 


y fdx: 2 


\£;V 


dxi 


(65) 


It allows us to use the following well-known equations and definitions. The metric tensors of 
the Euclidean space in curvilinear coordinates x k are 


9ki{x) = Hi Ski, g k \x ) = Ski, 


1 


HI 


and indices can be raised and lowered by this metric 

u = g ( x)ui , u k = g M u . 

For an orthogonal basis, we have 


n=3 n=3 


9 = 


det(g k i) = \\g kk = Y1 H ^ 


and 


k=1 k=1 


n =3 


J =y m\ = UHk = H 1 H 2 H 9 . 

k=1 

Then the volume differential 3-form is given by 

vol Q = \/\g{x)\ dx 1 A dx 2 A dx 3 . 

The correspondent volume for arbitrary set W in local coordinates is 


vol a (W) = / \/\g{x)\ dx 1 dx 2 dx 3 . 

Jw 


( 66 ) 


(67) 


( 68 ) 


(69) 


(70) 


(71) 
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Using the metric tensor () 66 lh the gradient of scalar function p is the vector held 


n =3 


grad (Q) (p= ^e k g 


k,l =1 


kir r) &P_ 


( 72 ) 


and the divergence of vector held u = u k e k is defined as the scalar function by 


n =3 


div (a) U = ]T 


d{\f\g{x)\u k ) 


^ VW)\ dxk 


(73) 


The Laplace-Bcltrami operator is defined as the divergence of the gradient. Using the definitions 
of the gradient and divergence, the Laplace-Beltrami operator applied to a scalar function p is 
given in local coordinates by 


n =3 


A {a) (p = div (a) grad (a) <p = '*T, 


1 _ d_ 

ut VW\ dxk 


V\g(x)\g kl {x)^ 


dx l 


Using the effective coordinates X k = X k (a k ,x k ), which are dehned in 
the nabla operator (the del operator) by 


d 


£ 


d 


dX k ci(a k ,x k ) dxk 


(74) 

we can define 

(75) 


where (a) = ( 01 , 0 ( 2 , 0 ( 3 ) is the multi-index, and ci(a k ,x k ) is dehned by fjMjl . We can use 
the well-known relations for the gradient, divergence, the curl operator, the scalar and vector 
Laplace operators through the Lame coefficients. The gradient for scalar held 


fc=i 


grad (Q) p = 

The divergence for vector held 

3 

div( Q ) u = 


1 dp 
Hi dx k 




(76) 


d / H, H , H, 


k= 1 


H l H 2 H,dx k 


H h 




(77) 


The curl operator for vector held 


curl (a) u = ^ 


k,i,j =1 


//, h 2 h, 


' ^Ijk H k 


d{H k u k ) 
dxj 


(78) 
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where e^ k is the Levi-Civita symbol. The scalar Laplacian 




H\ H 2 II, dx k 


d f Hi H, /hi Op 
Hi dx k 


(79) 


fc=1 - - -—" ' 

Using that H k = ci(«fc,Xfc), we get dH k /dxi = 0 for k ^ l, and the divergence 


div (a) u = ^ 


k=1 


1 du k 
Ilk dx k ' 


The curl operator for vector field 


r-r 1 du k 
curl (a) u = 2^ jf e i e ijk 


H, 

k,i,j=1 J 


dxj’ 


where e^k is the Levi-Civita symbol. The scalar Laplacian 


(80) 


(81) 


A <*)^ = 52 


1 d f 1 dp 


k=1 


Ilk dx k V Hk dx k 


(82) 


by 


We can define a differential operator that takes into account the density of states ci(a k) x k ) 


d 


d 


d 




dX k ci (a k pxk) dx^ 

where Ci(cxfc,a;fc) corresponds to the non-integer dimensionality along the W-axis and it is 
dehned by (1641) . These derivatives cannot be considered as derivatives of non-integer orders (as 
fractional derivatives) or as fractal derivatives (derivatives on fractal set). The operators d Xk>ak 
are usual differential operators of the first order that is dehned for differentiable functions on 


The form of derivatives (1551) is analogous to differential operator suggested in [39]. The 
main difference is that we define the operators (1531) for the non-integer dimensional spaces in 
the framework of the measure product approach that is described by Palmer and Stavrinou 
[67 ] (see also [39] SHUT]). The differential operators suggested in [39] are dehned for modified 
Riemann-Liouville fractional integral of orders a k in the integer dimensional space. 

Using the operators ([53]) . we can introduce generalized vector differential operations. The 
gradient 

3 

grad (Q) <p(r) = ^ e fc d Xkt0lk p(r), (84) 

k= 1 

where e k are unit base vector of the Cartesian coordinate system. The divergence of the vector 
held u(r) = e k u k ( r) is 

3 

div (a) u(r) = d XkiClk u k ( r). (85) 

k=1 
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( 86 ) 


The curl for the vector field u(r) = e k u k ( r) is 

3 

curl (a) u(r) = ^2 e i£ikid Xk , ak Ui( r), 
k,i,l=1 


where e^i is the Levi-Civita symbol (or alternating symbol). 

Using (184)) and (1851) we can define the second order differential operators such as the scalar 
Laplacian and vector Laplacian. The scalar Laplacian has the from 

5 A (Q) (^(r) = div (a) grad (Q) r). (87) 

The vector Laplacian |l04j has the from 

V A (a) u(r) = grad (o) div (a) u(r) - curl (a) curl (a) u(r). (88) 

The scalar Laplacian (187[) can be represented by using the usual partial derivatives 


’ A( a) 9?(r) = 


3 


k 


—[ c\(oi k ,x k ) 


d 2 if 

dxl 


a k - 1 d(p \ 
x k dx k J ’ 


(89) 


where Ci(a k ,x k ) is density of states (164]) along the X fc -axis for model with the non-integer 
dimensional space. It is easy to see that this operator does not coincide with the Laplace 
operators proposed by Palmer and Stavrinou [67], by Stillinger [61], and by Calcagni [54] • The 
main advantage of the suggested Laplace operator (j89|) is that firstly, it is obtained as the action 
of the gradient and divergence, and secondly, it is adapted for models with non-integer spatial 
dimensions. 

The differential operators of the first order such as the gradient (1761) . the divergence (1771) . 
the curl operator f[78|l . and the second order differential operators such as the scalar Laplacian 
f[87p and the vector Laplacian (ISSlh allow us to describe anisotropic fractal media and materials 
in the framework of continuum models with non-integer spatial dimensions. 


5 Examples of application 


5.1 Poisson’s equation 


Let us consider the Poisson’s equation for a fractal medium that is distributed along the positive 
half-X-axis 

S ^{a)f(x) = f(x). (90) 

Here we use the Laplace operator suggested in this paper. Equation (j90l) for single-variable 
case can be written as 


1 / d 2 (p(x) 

cl(a,x) \ dx 2 



f(x) (x > 0). 


(91) 
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The general solution of this equation is 


<p{x) = Ci + C 2 x a - a(r( ^ /2 ))2 (/ ^ 1 dx ~ x ° J f( x ) x ° 1 dx ) ■ ( 92 ) 

Let us compare this result with the solution of Poisson’s equation Laplace operator /C 2 . The 
Poisson’s equation of the form 

= f(x) (93) 

has the general solution 

<p(x) = C l X ^- a) / 2 + C 2 X {l - a) / 2 +x^ /2 J f(x)x {1 ~ a)/2 dx-x {1 - a)/2 J f{x)x (1+a)/2 dx. (94) 
In equation f[93lh we use the Laplace operator (I89|) for single-variable case. 


5.2 Timoshenko beam equations for fractal material 

The Euler-Bernoulli beam theory gives a simplification of the linear theory of elasticity, which 
provides a means of calculating the load-carrying and deflection characteristics of beams and 
it covers the case for small deflections of a beam, which is subjected to lateral loads only. Note 
that the Timoshenko beam equation for fractal materials is discussed in the papers [34, 39]. 

In the Timoshenko beam theory without axial effects, the displacement vector u(a;, y, z, t ) 
of the beam are assumed to be given by 


u x (x,y,z,t) = -zip(x,t) u y (x,y,z,t) = 0, u z (x,y,t) = w(x,t), 


(95) 


where (x, y, z ) are the coordinates of a point in the beam, u x , u y , u~ are the components of the 
displacement vector u , (p = <p(x,t) is the angle of rotation of the normal to the mid-surface of 
the beam, and w = w(x,t) is the displacement of the mid-surface in the ^-direction. 

Using a model with non-integer dimensional space for fractal media, we can use the deriva¬ 
tives 


&x,a 


Ci 1 


(a x,x)Dl, d"=(d„) n (neN), 


(96) 


where Ci(a x ,x) is defined for non-integer dimensional space by 


Cl (ci!x, •^) 


Tj-aic/2 

rp) 


rQ I &X 1 


(97) 


instead of the usual derivatives D l x and D™ for fractal materials. If we use the derivatives (j96]h 
then the Timoshenko equation for fractal beam can be derived from the force and moment 
balance equations 


p A = a,„Q. p /"‘ l Dpf = Q- d„M, 


(98) 


with the bending moment 

M = -EI^d x , a ip, 


(99) 
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and the shear force 


( 100 ) 


Q — kG A (d X:a w - pj 


Here is the second moment of the fractal beam’s cross-section. 

The Timoshenko equations for homogeneous fractal beam has the form 


p A D 2 w = k G Ad Xj0l (d Xta w — p), (101) 

pI& AV = kGA (d x , a w - ip) + E J (d) d 2 x>a p. (102) 

The Timoshenko fractal beam equations (1101 j) and (11021) also can be derived from the 
variational principle. The Lagrangian of the Timoshenko fractal beam has the form 

Cgtfb = \p lid) (D\ip{x,t )) 2 + \pA ( D\w(x,t )) 2 - 

- ^( kGA ) (d Xta w(x, t ) - (p(x, t)) 2 - ^( EI {d) ) (d Xta (p(x, t)f . (103) 

The stationary action principle gives the equations 


— -D 1 

dw t 





0, 


(104) 


dC 

dip 


~D] 


DC 


dD 1 


iv 


Dl 


dC 

dD^ 


= 0. 


(105) 


Equations (11041 - 11051) are the Euler-Lagrange equation for the model of fractal material de¬ 
scribed by the Lagrangian (11031) . Substitution of (I103|) into equations (11041 - 1105(1 gives the 
Timoshenko fractal beam equations (11011) and (11021) that can be presented as 


p A D'fw 


kGA d x . Q 



w 




(106) 


pl {d) D 2 t ip = kGA (d x>a w - p) +EI& dl a p. (107) 

If a = 1 then equations (I106H107]) are the Timoshenko equations for beam with homogeneous 
non-fractal material. 

For models with non-integer dimensional spaces, solutions of equations for fractal materials 
can be obtained from solutions of equations for non-fractal materials. Let w c (x,t) and ip c (x,t) 
be solutions of (1106111071) with a = 1 and x > 0, i.e., the Timoshenko equations for homogeneous 
non-fractal beam. Then solutions wp(x,t) and (pp(x,t) of equations (jl06H107[) for fractal beam 
with 0 < a < 1 can be represented by 

”' F(l ' f) = ro '(f^j W “’' 1,i )’ Mx ’* ) = ^(r ' (108) 
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5.3 Euler-Bernoulli fractal beam 


As an example, we consider the equation for the Euler-Bernoulli homogeneous fractal beam in 
the absence of a transverse load (q(x) = 0), 


pAD 2 w(x,t) + E J (d) d^ a w(x,t) = 0. 


(109) 


This equation can be solved using the Fourier decomposition of the displacement into the sum 
of harmonic vibrations of the form w(x,t) = Re[w(x) exp(— iut)], where u is the frequency of 
vibration. Then, equation (11091) gives the ordinary differential equation 


— pAu 2 w(x ) + E1^ d^ a w(x) = 0 . 

The boundary conditions for fractal beam of length L fixed at x = 0 are 

w(0) = 0, (d* ia w)(0) = 0, 

( dl,a W )( L ) = 0 , ( d la W )( L ) = 0 - 

The solution for the Euler-Bernoulli homogeneous fractal beam is defined by 
WF,n{x ) = Wo (cosh(fc n a; a ) — cos (k n x a ) + C n (a) (sin (k n x a ) — sinh(/c n a; a ) 


where Wo is a constant, and 

cos (k n L a ) + cosh [k n L c 


( 110 ) 

( 111 ) 

( 112 ) 


x G [0; L\, 

(113) 


CJa) = 


fort. — 


7r 


a/2 


P Auj l \ 


2 \ 1/4 


sm(k n L a ) + sinh(fc ?l L a ) ’ ” T(a/2) \ E E d ) J 

For boundary conditions (111111112]) . the solution (II 13j) exist only if k n are defined by 

cosh (k n L) cos (k n L) + 1 = 0. 


(114) 


(115) 


This trigonometric equation is solved numerically. The corresponding natural frequencies of 
vibration are u n = k 2 \J(E /( d ))/p A. For a non-trivial value of the displacement, w o has to 
remain arbitrary, and the magnitude of the displacement is unknown for free vibrations. Usually 
Wq = 1 is used when plotting mode shapes. 


6 Conclusion 

We give a review of possible approaches to describe anisotropic fractal media. We focused on 
two approaches based on the fractional space and the non-integer dimensional space. There 
approaches allow us to describe anisotropic fractal media and materials in the framework of con¬ 
tinuum models by using the concept of density of states [27] and the product measure method. 
Fractal medium is considered as a medium with non-integer mass dimension. The non-integer 
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dimensionality is a main characteristic property of fractal materials. Therefore we suggest an 
application of differentiation and integration over non-integer dimensional spaces as natural 
way to describe fractal media [681 |69j. Although, the non-integer dimension does not reflect 
all specific properties of real fractal materials, it allows us to formulate continuum models to 
derive important conclusions about the behavior of the media. In this paper a generalization of 
the vector calculus for multi-fractional and non-integer dimensional spaces is proposed as tools 
to describe anisotropic fractal media and materials in the framework of continuum models. We 
suggest a generalization of vector calculus for non-integer dimensional space that is product of 
spaces with different dimensions. The product measure method allows us to describe anisotropic 
fractal materials by taking into account various non-integer dimensions in different directions. 
The differential operators of the first order such as the gradient (1761) . the divergence (1771) . the 
curl operator (1751) . and the second order differential operators such as the scalar Laplacian (157|) 
and the vector Laplacian (1881) . are suggested to describe anisotropic fractal media and mate¬ 
rials by continuum models with non-integer dimensional spaces. To demonstrate some simple 
applications of proposed approach to the description of fractal materials, we consider the Pois¬ 
son’s equation for fractal medium, the Euler-Bernoulli fractal beam and the Timoshenko beam 
equations for fractal material. 
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